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Derealization Transition and Localization Length 
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One dimensional system of Dirac fermions with a random- varying mass is studied by the transfer- 
matrix methods which we developed recently. We investigate the effects of nonlocal correlation of 
the spatial- varying Dirac mass on the derealization transition. Especially we numerically calculate 
both the "typical" and "mean" localization lengths as a function of energy and the correlation length 
of the random mass. To this end we introduce an imaginary vector potential as suggested by Hatano 
and Nelson and solve the eigenvalue problem. Numerical calculations are in good agreement with 
the results of the analytical calculations. 



PACS: 72.15.Rn, 73.20.Jc, 72.10.Bg 

One of the most important problems in condensed mat- 
ter physics is the localization phenomena in random- 
disordered systems 0. At present it is believed that 
all states tend to localize in disorder systems in two 
and lower dimensions. In some special cases, however, 
some specific states remain extended even in the pres- 
ence of strong disorders. System of Dirac fermions with 
a random- varying mass in one dimension has been stud- 
ied from this point of view |^ |^]. In the previous papers 
P-|Tl| we studied the effect of nonlocal correlation of the 
random mass on the extended states which exist near 
the band center. For numerical studies, we reformulate 
the system by transfer-matrix formalism, and obtained 
eigenvalues and wave functions for various configurations 
of random telegraphic mass . We verified that the den- 
sity of states obtained by the transfer-matrix methods is 
in good agreement with the analytical calculation in Ref. 
p0| . In this paper we shall introduce an imaginary vec- 
tor potential into the system of the random-mass Dirac 
fermions and study a localization-delocalization phase 
transition by varying the magnitude of the vector po- 
tential. Through this study we shall obtain the localiza- 
tion length of the states as a function of energy and the 
correlation length of the random mass. This method of 
calculating the localization length is based on the idea 
by Hatano and Nelson |j[2| . 

We shall consider a Dirac fermion in one spatial di- 
mension with a coordinate-dependent mass m(x) and in 
an imaginary vector potential g, whose Hamiltonian is 
given by, 



Tt = j dxip^hip, 
h = -icr z (d x +g) + m(x)(T y , 



(1) 
(2) 



where a are the Pauli matrices. This fermion model is 
an low-energy effective model of random-hopping tight 
binding models and random-bond spin chains |J. We 
introduce the components of ip as tp = (u,v). In terms 



of them the Dirac equation is given as, 

d_ 

dx 



+ g+ m(x) I u(x) = Ev(x), 



-— — g + m{x) ] v{x) = Eu(x). 



(We follow the notations in Ref. From Eqs. 
obtain the Schrodinger equations, 



d 2 d 

"fai ~ 2 9 fa ~ m '( x ) + ( m ( x ) - g 2 ) ) u i x ) 



(3) 



we 



E 2 u(x), 



(4) 



and similarly for v(x). 

In this paper we restrict the shapes of m(x) to multi- 
soliton-antisoliton configurations |H| . The multi-soliton- 
antisoliton configurations are given by, 



1) 



i{x) — m(Q(x — oti) 

i 

+ ^m(0(-x + /%)-l), 



(5) 



where ai's(/3j's) are positions of solitons(anti-solitons). 
An example of m(x) is given in Fig.l. If we vary the 
distances I between soliton and anti-soliton according to 
the exponential distribution like, 



P(l) =Xexp(--L) 
2A V 2A/ 



(6) 



where A is a parameter, then m(x) has the following cor- 
relation H], 

A 

[ m(x) m(y) ] ons = ^ ex P (-\ x - J/IA), (7) 
A 



where y y essentially corresponds to the height of the 
soliton and anti-soliton, i.e., fh in (||). From (^), A is 



1 



the correlation length of the random mass and the limit 
A — > corresponds to the white-noise case. In subse- 
quent papers, we shall study Dirac fermions with long- 
range correlated random mass by using the methods ex- 
amined in this paper Q. There we expect some interest- 
ing phenomena like existence of nontrivial mobility edge, 
nonuniversality of the multi-fractal exponents, etc. Stud- 
ies in this paper show that the imaginary- vector-potential 
methods for calculating the localization lengths are reli- 
able and we shall use them for studies on random systems 
with long-range correlated disorder. 



FIG. 1. An example of configurations of solitons and 
anti-solitons. 



For the vanishing imaginary vector potential, we 
solved the Schrodinger equations in (Q) under the pe- 
riodic boundary condition with various multi-soliton- 
antisoliton configurations of m(x) using the transfer- 
matrix method, and obtained the energy spectrum and 
wave functions ||. 

Effect of the imaginary vector potential was discussed 



by Hatano and Nelson |12|. Let us denote the eigenfunc- 
tion of energy E for g = as ^q(x), and suppose the 
shape of ^0(2;) as 



V (x) = cxp 



Co 



(8) 



where £0 is the localization length and x c is the center of 
this localized state. When we turn on the constant imag- 
inary vector potential g, the eigenfunction is obtained 
from ^o(x) by the "imaginary" gauge transformation, 



1 I C\ / \ 

w(z) ^ exp g (x - Xc) 



(9) 



This means that the localization length of this eigenfunc- 
tion is 



C 9 = 



Co 



1 + <?Co 
Co 

i-<?a> 



(X > X c ) 
(X < X c ). 



(10) 



At the point g = 1/Co: localization length for x < x c 
diverges, and if the imaginary vector potential g is in- 
creased more, £ g for x < x c becomes negative, and this 
eigenfunction cannot satisfy the periodic boundary con- 
dition even if the length of system is large enough. So if 
localized eigenstate with energy E disappears at g = g c 
as g is increased, then the localization length of the eigen- 
state ^o(x) is l/<7c- Actually it is shown that for g > g c , 
energy eigenvalue of the state has an imaginary part and 
the state is extended as we shall see shortly [Q . 

The transfer-matrix methods can be easily extended 
for the case of nonvanishing g. We obtain energy eigen- 
values and cigenfunctions numerically. First of all, we 
show that the derealization transition actually occurs. 
In Fig. 2, we show the wave functions of a low-lying state 
in vanishing and nonvanishing imaginary vector poten- 
tial. For g = 0, the state is obviously localized whereas 
at g = 0.03 the state becomes extended and the en- 
ergy eigenvalue has an imaginary part. As discussed in 
Rcf. [^2|, the density distribution of a particle is given 
by \^>(x, — g)^(x, g)\ where ^{x, —g) is equal to the left 
eigenfunction. 
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FIG. 2. An example of the localized and the extended wave 
functions in vanishing and nonvanishing imaginary vector po- 
tential g = and g = 0.03. The state of g = 0.03 has a 
complex energy eigenvalue. 



Let us turn to the localization length. For the white- 
noise case [m(x)m(y)] ens = A S(x — y), "typical" local- 



2 



ization length or the inverse of the Lyapunov exponent 
was obtained as pJl5f 



&(E) = \\nE/2A\ 



(11) 



Numerically the typical localization length is ob- 

tained by averaging over localization lengths of all eigen- 
states with energy E. On the other hand, Balents and 
Fisher calculated the averaged Green function and ob- 
tained the mean localization length from the spatial de- 
cay of the Green function j5| . The result is 



U(E) = \\nE/2A\ 



(12) 



Case of nonlocally correlated random mass was studied 
in Refs. |Io| , |lT|] and €, m {E) is obtained as a function of A 
inEq.(@). 

By numerical calculation we obtain both the typical 
and mean localization lengths. As we mentioned above, 
the typical localization length is the average over all solu- 
tions of the Schrodinger equation (Q), whereas the mean 
localization length is determined by the states which 
make dominant contributions to the Green function, i.e., 
which have large localization length. 

The result of the numerical calculation of the typical lo- 
calization length of the white-noise case is given in Fig. 3. 
We show the ratio of the numerical results to the analyt- 
ical calculation in Eq. ([ll]) in order to compare these two 
results. Therefore if the energy dependences of the lo- 
calization lengths obtained numerically and analytically 
are the same, this ratio should be constant. In Fig. 3, the 
ratio seems constant over the whole range of E. 
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FIG. 3. The comparison between the analytical and numer- 
ical results of the typical localization length for the (almost) 
white-noise case, i.e., small-A case. In the numerical calcula- 
tion, we set L(system size) = 50, A = 1/60, and the energy 
slice 8E = 0.02. This result is averaged over 2500 trials. The 
analytical and numerical results are in good agreement. 



In Fig. 4, we show the numerical results of the "mean" 
localization length. Here we use the solutions to the 
Schrodinger equation which have long localization length. 
More precisely, "large" localization length £ means 



the one which satisfies £ > (the "typical" localization 
length)+ 1.5 a in each energy slice. From Fig. 4, we can 
conclude that the energy dependence of the mean local- 
ization length obtained numerically is in agreement with 
Eq@. 
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FIG. 4. The comparison between the analytical and numer- 
ical results of the mean localization length. We use the data 
of the states which have large localization length. (See the 
text.) We show the ratios of the numerical calculations both 
to Eqs.(|ll"|) and (|l^). The ratios are normalized at E = 0.19. 



From the above studies, we can conclude that the 
above methods of calculating the localization lengths are 
reliable. 

We shall turn to the case of the nonlocally-correlated 
disorder. In the white-noise limit (A = case in Eq.(^)), 
the localization lengths diverge only at E = 0, that is, 
extended states exist only at E = 0. If we let A > 0, 
the random mass becomes nonlocally-correlated, and the 
critical energy or the mobility edge at which the dereal- 
ization transition occurs may change. 

We investigate the "typical" and "mean" localization 
lengths in the case of nonvanishing A's. The behaviour of 
the "typical" and "mean" localization lengths obtained as 
in the white-noise case are given in Figs. 5 and 6. It seems 
that there is no A-dependence in the typical localization 
length. On the other hand, Fig. 6 shows that the mean 
localization length has a small but finite dependence on 
A. 
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FIG. 5. The behaviour of the localization length in the case 
of nonvanishing A's. The ratio is the "typical" localization 
length to | laE\. The ratio is normalized at E — 0.19. 
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FIG. 6. The behaviour of the localization length in the 
case of nonvanishing A's. The ratio is the "mean" localization 
length to | In | 2 . The energy slice 8E = 0.03, and the ratio 
is normalized at E — 0.19. 



From the above calculations, we conclude that the ef- 
fect of the short-range correlations in disorders is not so 
large. Especially the result indicates that the derealiza- 
tion transition occurs at E < 0.03. (If the mobility edge 
exists at E c > 0, the ratio in Fig. 5 or 6 must diverge 
at E c .) The derealization transition probably occurs at 
E = 0. 

In the previous paper |n| we calculated the localiza- 
tion length for the random mass with the short-range cor- 
relation (^) . We obtained the "mean" localization length 
to the 1st order of A by means of the Green function 
method. The result is 



1 /In 

m = - A 



|_E |2 
\1A\ 



- 41 In — 



+ 0(A 2 ). (13) 



In Fig. 7 we show the ratios of the numerical result to the 
analytical calculation up to the 0th and the 1st order of A. 
This shows that the analytical result with the 1st order 
correction of A is in better agreement with the numerical 
result, but the correction is small. 
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FIG. 7. The comparison between the numerical result and 
the analytical one in the case of relatively large A. We show 
the ratio of the localization length calculated numerically to 
the ones obtained analytically. Here we used the analytical 
results of the 0th and the 1st order of A. Here we set A = 1/4 
and L = 50. The energy slice SE = 0.03, and the ratio is 
normalized at E = 0.19. 



From the investigations given as far we can conclude 
that the numerical methods used in this paper are reli- 
able for calculating the localization lengths. It is very 
interesting to study the case of disorders with a long- 
range correlation. We expect that a nontrivial mobility 
edge E c > exists for a certain long-range correlated 
random mass. Actually the one-dimensional Anderson 
model with long-range correlated disorder was studied 
|l6| , and it is shown that there exists a nontrivial mobil- 
ity edge. 

We can also calculate exponents of the multi-fractal 
scaling [17J by the transfer- matrix methods ||. These 
problems are under study and results will be reported in 
a future publication |14[ . 
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